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ABSTRACT 


Earthquake  data  from  Central  California  over  the  period 
January  1,  1969,  to  December  31,  1971,  and  from  Western 
Nevada  over  the  period  August  31,  1954,  to  December  31,  1961, 
have  been  analyzed  to  investigate  the  hypothesis  that  the 
Earth  tides  may  be  the  triggering  mechanism  for  the  occurrence 
of  earthquakes.   One  method  used  compared  the  tidal  components 
and  the  component  rates -of -change  at  the  time  and  location  of 
actual  earthquakes  to  those  for  random  events  uniformly 
generated  over  the  same  time  period.   A  second  method  com- 
pared the  differences  in  the  Earth  tides  computed  across  the 
breadth  of  a  fault  zone  at  the  time  of  actual  earthquakes  to 
the  differences  at  the  times  of  random  events  uniformly 
generated  over  the  same  time  period.   Statistical  tests  were 
performed  to  see  if  the  distributions  from  the  actual  events 
were  the  same  as  the  corresponding  distributions  from  the 
randomly  generated  events. 
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I.   INTRODUCTION 

Every  point  on  the  surface  of  the  Earth  is  subject  to 
three  types  of  forces :   the  force  of  gravity  due  to  the 
attraction  of  the  mass  of  the  Earth,  the  centrifugal  force 
due  to  the  rotation  of  the  Earth,  and  a  force  due  to  the 
attraction  of  extra-terrestrial  bodies,  mainly  the  Sun  and 
the  Moon.   The  resultant  forces  on  the  surface  of  the  Earth 
are  generally  referred  to  as  the  Earth  tide  forces. 

Because  of  the  stresses  that  these  Earth  tides  create  in 
the  crust  of  the  Earth,  it  has  been  long  felt  that  they  may 
be  a  significant  factor  in  the  triggering  of  earthquakes. 
Several  previous  studies  have  been  undertaken  in  which 
attempts  were  made  to  correlate  the  gravitational  component 
of  the  Earth  tides  with  the  occurrence  of  earthquakes.   In 
one  study,  Knopoff  [Ref.  1]  concluded  that  there  was  no 
correlation  between  this  tidal  component  and  the  occurrence 
of  the  3800  earthquakes  in  Southern  California  that  he  inves- 
tigated.  Shlien  [Ref.  2]  came  to  a  similar  conclusion  in  his 
study  of  several  widely  separated  regions  of  the  world. 

A  team  of  Russian  researchers  [Ref.  3]  contend  that  these 
and  other  similar  studies  have  suffered  from  two  serious 
drawbacks.   First,  they  considered  only  the  gravitational 
component  of  the  Earth  tides.   More  important  may  be  the 
tangential  components  which  depend  very  strongly  on  the 
orientation  of  a  fault  in  a  region.   In  general,  the  extremes 
in  the  tangential  components  do  not  occur  at  the  times  of  the 
extremes  of  the  gravitational  component.   Second,  they  tried 


to  cover  too  large  an  area.   Knopoff  considered  a  region  of 
18  square  degrees,  and  the  regions  in  Shlien's  study  varied 
between  50  and  450  square  degrees.   Since  tides  are  time 
dependent  and  periodic  in  nature,   considerations  of  large 
regions  have  averaging  effects  on  their  variation,  possibly 
masking  dependencies  present. 


II.   OBJECTIVE 

The  objective  of  this  study  was  three-fold.   Given  a 
region  of  investigation  and  a  specific  time  period:   (1) 
compare  the  tidal  forces  from  a  series  of  actual  earthquakes 
to  these  forces  from  a  series  of  random  events  generated 
uniformly  over  the  same  time  period,   (2)  compare  the  rates- 
of-change  of  the  tidal  forces  from  the  actual  events  to  those 
from  the  randomly  generated  events,  and  (3)  compare  the 
differences  in  each  of  the  tidal  forces  over  the  breadth  of 
the  region  from  the  actual  events  to  those  from  the  randomly 
generated  events.   These  comparisons  were  used  to  test  the 
hypothesis  that  the  distribution  of  a  force,  rate,  or  differ- 
ence from  the  actual  events  was  the  same  as  the  corresponding 
distributions  from  the  randomly  generated  events.   Rejection 
of  the  hypothesis  would  indicate  a  difference  in  the  distri- 
butions being  compared.   Thus,  an  evaluation  can  be  made  as 
to  whether  a  significant  dependency  exists  between  a  force, 
rate,  or  difference  and  the  occurrence  of  earthquakes  for  the 
particular  region  being  studied. 


III.   ORGANIZATION 

For  the  purpose  of  this  analysis,  the  tidal  forces  were 
resolved  into  five  quantities.   The  orientation  and  con- 
struction of  these  quantities  were  as  follows:   The  gravita- 
tional component  was  formed  along  the  radius  of  the  Earth  at 
the  earthquake  epicenter,  the  normal  component  was  formed 
perpendicular  to  both  the  gravitational  component  and  the 
strike  of  a  fault,  the  axial  component  was  formed  perpendic- 
ular to  the  gravitational  component  and  parallel  to  the 
strike  of  a  fault,  the  tangent  magnitude  was  the  resultant 
formed  by  the  combination  of  the  magnitudes  of  the  normal 
and  axial  components,  and  the  total  magnitude  was  the 
resultant  formed  by  the  combination  of  the  magnitudes  of  the 
gravitational,  normal,  and  axial  components.   By  observing 
how  these  five  quantities  changed  over  a  small  change  in 
time,  the  rate-of -change  of  each  quantity  also  was  approxi- 
mated.  By  calculating  these  tidal  quantities  on  both  sides 
of  a  fault  zone,  the  differences  in  these  quantities  also 
were  calculated. 

It  should  be  noted  here  that  the  foregoing  definitions  of 
the  normal  and  axial  components  of  the  tidal  forces  usually 
precludes  their  being  exactly  tangent  to  the  surface  of  the 
Earth  owing  to  the  fact  that  the  Earth  in  not  a  perfect 
sphere.   Since  this  angular  deviation  from  tangency  was  so 
small  and  was  nearly  constant  over  the  small  range  in  lati- 
tude of  the  regions  investigated,  no  correction  was  made  in 
calculating  component  values.   This  correction  should  cancel 
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out  during  the  comparisons  of  component  frequency 
distributions  and,  therefore,  have  little  influence  on  the 
results  of  the  comparisons. 

In  selecting  regions  for  investigation,  zones  of  high 
seismic  activity  were  desired  so  that  relatively  small  regions 
could  be  selected  and  still  yield  a  sufficient  number  of 
earthquake  occurrences  to  support  a  detailed  analysis.   The 
regions  chosen  were  small  sections  of  the  San  Andreas  fault 
in  Central  California  and  the  Dixie  Valley  fault  in  Western 
Nevada,  both  of  which  are  located  in  the  Western  United 
States  (Figure  1) .   These  two  regions  were  chosen  because  of 
their  differing  fault  structure.   The  San  Andreas  fault  is  a 
strike-slip  fault  where  the  land  mass  to  the  southwest  side 
of  the  fault  moves  northward  in  relation  to  the  land  mass  on 
the  northeast  side  of  the  fault  (Figure  2).   In  the  partic- 
ular region  of  the  San  Andreas  fault  investigated,  the  strike 
of  the  fault  is  approximately  40  degrees  west  of  north  as 
measured  from  longitude  121  W.   The  Dixie  Valley  fault  is  a 
graben  fault  where  the  material  in  the  center  of  the  fault 
zone  moves  downward  in  relation  to  the  sides  of  the  fault 

(Figure  3) .   The  strike  of  this  fault  is  approximately 
parallel  to  longitude  118  W.   Co-ordinates  of  these  regions 
are  given  in  Table  1. 

In  choosing  a  time  period  over  which  earthquakes  were  to 
be  analyzed,  a  large  enough  time  period  was  needed  so  that 
bias  due  to  the  periodic  fluctuation  of  the  Earth  tides 

(monthly,  annually,  etc.)  would  be  damped  or  averaged  out. 


This  was  needed  so  that  differences  found  in  the  comparison 
of  these  quantities  would  not  be  attributable  to  the  cyclic 
nature  of  the  tides.   For  the  specified  portion  of  the  San 
Andreas  fault,  1732  earthquakes  of  Richter  magnitude  0.0  or 
greater  occurred  in  the  three  year  period  from  January  1, 
1969,  to  December  31,  1971.   For  the  portion  of  the  Dixie 
Valley  fault,  the  available  data  prior  to  1962  was  for  shocks 
of  magnitude  4.0  or  greater,  and  from  1962  on,  for  shocks  of 
magnitude  0.0  or  greater.   So  as  not  to  unduly  bias  the 
distribution  of  data  with  the  increased  frequency  of  recorded 
small  shocks  after  1962,  only  the  time  period  through  1961 
was  considered.   During  the  seven  year  period  from  August  31, 
1954  (the  start  of  detailed  recording  in  this  area) ,  to 
December  31,  1961,  100  earthquakes  of  magnitude  4.0  or  greater 
occurred . 
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Strike-Slip  Fault  Structure 

Figure  2 


I 


Dixie  Valley 
Graben  Fault  Structure 

Figure  3 


12 


San  Andreas 


Dixie  Valley 


Corners  of  Area 


Center  of  Area 
Difference  Points 

N-S  Comparison 

W-E  Comparison 


121  W  18 

121  W  30 

120  W  48 

120  W  36 

121  W  03 

121  W  36 

120  W  36 

121  W  18 

120  W  48 

121  W  30 
120  W  36 


36  N  54 

36  N  42 

36  N  00 

36  N  12 

36  N  27 

35  N  54 

36  N  54 

36  N  54 

36  N  00 

36  N  42 

36  N  12 


117  W  48 

118  W  18 
118  W  18 

117  W  48 

118  W  03 

118  W  30 

117  W  30 

118  W  03 
118  W  03 

118  W  18 

117  W  48 


39  N  48 

39  N  48 

38  N  54 

38  N  54 

39  N  21 

39  N  21 

39  N  21 

39  N  48 

38  N  54 

39  N  21 
39  N  21 


Area  Co-ordinates 
Table  1 


13 


IV.   EARTH  TIDE  DISCUSSION 

The  actual  tidal  data  were  generated  through  the  use  of  a 
computer  routine  (unpublished)  by  Professor  R.  II.  Shudde, 
Naval  Postgraduate  School.   Much  work  has  been  accomplished 
over  the  past  few  years  on  the  subject  of  theoretical  Earth 
tides.   Melchior  [Ref.  4]  has  one  of  the  better  discussions 
on  the  subject.   He  goes  into  great  detail  concerning  the 
computation  of  the  theoretical  Earth  tides  using  an  equi- 
potential  surface  argument.   His  procedure  determines  the 
locations  of  the  Sun  and  the  Moon  in  relation  to  the  Earth  at 
a  particular  time  and  establishes  the  potential  surfaces  for 
the  tides.   Then  by  specifying  the  location  of  a  point  on  the 
Earth  in  relation  to  these  surfaces,  he  determines  the  tides 
at  that  location.   Shudde' s  program  takes  a  different  approach 
It  also  determines  the  locations  of  the  Sun  and  the  Moon  using 
the  procedures  in  reference  5,  but  then  at  a  particular  point 
on  the  Earth,  the  tides  are  calculated  directly  using  the 
standard  gravitational  force  laws.   For  a  more  detailed 
description  of  the  Earth  tide  calculations,  see  the  discussion 
in  Appendix  A.   It  will  suffice  here  to  say  that  all  Shudde ' s 
program  requires  for  inputs  are  the  time,  date,  and  geographic 
co-ordinates  of  an  earthquake  occurrence.   For  this  analysis, 
the  times  and  locations  for  the  actual  occurrences  were  taken 
from  references  6-8  and  from  data  supplied  by  the  National 
Center  for  Earthquake  Research.   The  dates  for  the  random 
occurrences  were  generated  uniformly  over  the  time  periods 
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specified  in  Section  III.   The  location  for  these  random 
events  was  taken  to  be  the  center  of  the  region  being 
investigated. 
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V.   ANALYSIS 

The  objective  of  this  analysis  was  to  test  the  hypothesis, 

HQ:  A  =  R 
where  A  is  the  distribution  of  data  from  the  actual  events, 
and  R  is  the  distribution  of  the  corresponding  data  from  the 
randomly  generated  events.   Rejection  of  this  hypothesis  for 
a  particular  force,  rate,  or  difference  would  indicate  a 
difference  in  the  two  distributions  being  compared  and,  hence, 
a  dependency  between  that  quantity  and  the  occurrence  of 
earthquakes . 

The  primary  statistical  test  used  to  test  this  hypothesis 

was  the  Kolmogorov-Smirnov  (K-S) ,  Two-Sample  test.   This  is  a 

non-parametric  test  independent  of  distributional  assumptions 

about  the  data  other  than  the  distribution  be  continuous.   It 

is  an  exact  test  and  operates  by  the  comparison  with  the 

actual  data  rather  than  with  grouped  data.   A  second  test 

used  was  the  Chi-Square  Goodness -of -Fit  test  which  operates 

by  the  comparison  of  grouped  data. 

The  K-S  statistic  "possesses  the  advantage  that  it  is  an 
exact  method,  whereas  the  Chi-Square  method  requires  a 
fairly  large  sample  to  justify  the  approximations  that  are 
needed  to  apply  it"  [Hoel,  1971]. 

The  Chi-Square  test  was  used  because  this  is  the  common 

statistical  test  used  in  many  studies.   The  K-S  test  was  used 

because  it  is  a  more  powerful  test  and  was  more  appropriate 

to  the  data. 

A.   KOLMOGOROV-SMIRNOV,  TWO-SAMPLE  TEST 

This  test  was  implemented  through  the  use  of  the  IBM 
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computer  subroutine,  K0LM2  [Ref .  10] ,  with  a  modified  sort 
routine  [Ref.  11].   Using  these  routines,  for  each  of  the 
five  tidal,  rate,  and  difference  quantities,  the  distribution 
of  data  from  the  actual  earthquakes  in  a  region  was  compared 
to  the  distribution  for  that  quantity  calculated  from  the 
randomly  generated  events.   The  two  sets  of  data  for  each 
distribution  were  each  sorted  into  increasing  order  and  formed 
into  sample  distribution  functions.   The  maximum  absolute 
ordinate  difference  between  the  two  distributions  was  deter- 
mined.  A  sample  size  adjustment  factor  was  calculated  by 
taking  the  square  root  of  the  ratio  of  the  sample  size 
product  to  the  sample  size  sum.   The  K-S,  two-sample  statistic 
Z,  was  then  calculated  by  taking  the  product  of  the  maximum 
difference  and  the  adjustment  factor.   That  is, 


Z  =  ™|Am-Sm|«/1^ 


where  A(x)  is  the  distribution  function  from  the  actual  events, 
S(x)  is  the  distribution  function  from  the  random  events,  N, 
is  the  number  of  events  in  A(x),  and  N^  is  the  number  of 
events  in  S(x).   The  use  of  the  adjustment  factor  is  necessary 
in  the  K-S  test.   If  not  applied  to  the  statistic  as  shown 
above,  the  inverse  of  the  factor  would  need  to  be  applied  to 
the  critical  values  listed  in  Table  2  to  keep  the  test  statis- 
tic consistent  with  the  critical  values.   In  this  case,  appli- 
cation of  the  adjustment  factor  to  the  statistic  simplified 
the  use  of  the  computer  routines  and  enabled  the  critical 
values  to  remain  constant  for  a  given  level  of  significance. 
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B.   CHI -SQUARE  GOODNESS-OF-FIT  TEST 

For  each  tidal  quantity,  the  maximum  range  covering  both 
the  tidal  data  from  the  actual  earthquakes  and  the  data  from 
the  randomly  generated  events  were  divided  into  20  equal 
intervals  and  the  data  was  grouped  into  these  intervals . 
Adjacent  intervals  were  collected  together  until  the  expected 
frequency  in  the  collected  intervals  exceeded  three  [Ostle, 
1963].   The  Chi-Square  statistic  was  then  calculated  as 

yl_  g  (Oi  -El)* 

where  0^  is  the  frequency  in  interval,  i,  from  the  actual 
events,  E^  is  the  expected  frequency  in  interval,  i,  and  m ' 
is  the  number  of  collected  intervals.   Selected  critical 
values  for  the  Chi-Square  test  are  given  in  Table  3. 
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Level  of  Value  of  Z  So  Large  as  to 

Significance  Call  for  Rejection  of  the  Hn 

.10  1.22 

.05  1.36 

.025  1.48 

.01  1.63 

.005  1.73 

Adapted  from  Beyer,  p.  429. 


Critical  Values  of  the 
K-S,  Two-Sample,  Statistic,  Z 

Table  2 
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Levels  of  Significance 
.10     .05    .025     .01    .005 
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17.3 

19.7 
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24.7 
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18.5 
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26.7 

28.3 

19.8 

22.4 

24.7 

27.7 

29.8 

21.1 

23.7 

26.1 

29.1 

31.3 

22.3 

25.0 

27.5 

30.6 

32.8 

23.5 

26.3 

28.8 

32.0 

34.3 

24.8 

27.6 

30.2 

33.4 

35.7 

26.0 

28.9 

31.5 

34.8 

37.2 

27.2 

30.1 

32.9 

36.2 

38.6 

28.4 

31.4 

34.2 

37.6 

40.0 

Adapted  from  Beyer,  p.  294. 


Critical  Values  of  the 
Chi-Square  Statistic 

Table  3 
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VI .  RESULTS 

In  the  analysis  of  the  tidal  data,  some  bias  in  the 
distributions  of  the  data  was  expected  due  to  the  physical 
distance  between  the  origin  of  the  actual  earthquakes  and  the 
geometric  center  of  the  region  at  which  the  random  events 
were  generated.   Although  this  bias  was  felt  to  be  minimal 
due  to  the  small  size  of  the  regions  and  the  large  number  of 
events,  additional  checks  were  performed  to  test  the  size  of 
this  bias.   Random  events  were  generated  at  the  north,  south, 
west,  and  east  extremes  of  each  region  and  the  five  tidal 
and  five  rate  quantities  were  calculated  at  each  location  for 
each  occurrence.   The  ten  quantities  at  the  north  and  south 
extremes  were  then  compared  using  the  K-S  test,  as  were  the 
ten  quantities  from  the  west  and  east  extremes.   The  results 
are  shown  in  Tables  4  and  5.   The  conclusion  was  that  bias  in 
the  data  due  to  the  physical  separation  of  the  events  was  so 
slight  as  to  have  an  insignificant  influence  on  the  final 
analysis.   It  must  be  stressed  here  that  the  significance  of 
the  following  results  pertains  only  to  the  particular  region 
and  the  nature  and  orientation  of  the  particular  fault  being 
analyzed. 

A.   SAN  ANDREAS 

The  K-S,  two-sample  statistics  for  each  of  the  five 
tidal,  five  rate,  and  five  difference  quantities  acting  in 
the  San  Andreas  region  are  shown  in  Table  6  and  the  Chi- 
Square  statistics  for  these  quantities  are  shown  in  Table  7. 
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Referring  to  Tables  2  and  3  for  a  0.05  level  of  significance, 
the  normal  tidal  and  gravitational  rate  quantities  were 
statistically  significant  under  both  tests.   None  of  the 
tidal  differences  was  statistically  significant  at  this  level 
in  either  test. 

B.   DIXIE  VALLEY 

The  K-S,  two-sample  statistics  for  each  of  the  five  tidal, 
five  rate,  and  five  difference  quantities  acting  in  the  Dixie 
Valley  region  are  shown  in  Table  8  and  the  Chi-Square  statis- 
tics for  those  quantities  are  shown  in  Table  9.   Referring  to 
Tables  2  and  3  for  a  0.05  level  of  significance,  the  normal 
rate  quantity  was  statistically  significant  under  both  tests. 
None  of  the  tidal  difference  quantities  was  significant  at 
this  level  in  either  test. 

After  much  of  this  analysis  had  been  completed,  a 
similarity  between  the  methods  of  calculating  the  rate  and 
difference  quantities  was  observed.   The  rate  quantities  were 
calculated  by  computing  the  tidal  quantities  at  a  location, 
allowing  a  small  unit  of  time  to  elapse,  and  then  computing 
the  tidal  quantities  again.   The  change  in  the  tidal  quantity 
divided  by  the  change  in  time  was  used  as  its  rate  of  change. 
Because  of  the  rotation  of  the  Earth,  changes  in  time  corres- 
pond to  changes  in  location  at  a  fixed  latitude.   A  lapse  of 
four  minutes  is  equivalent  to  a  distance  change  of  one  degree 
in  longitude.   For  small  time  lapses,  the  distance  considered 
moved  is  negligible.   For  larger  time  lapses,  the  calculations 
can  be  considered  as  being  at  two  separated  locations.  The 
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difference  quantities  were  calculated  by  computing  simultan- 
eously the  tidal  quantities  at  two  locations  separated  by  one 
degree  in  longitude.   Since  both  procedures  involve  the 
calculation  of  components  at  separated  locations,  there 
should  have  been  a  correlation  between  the  test  results  on 
the  rate  quantities  and  the  test  tesults  on  the  difference 
quantities.   Comparison  of  the  San  Andreas  results  in  Table 
6  showed  some  similarity  but  the  interpretation  was  clouded 
in  that  the  tangential  components  in  this  region  were  rotated 
through  an  angle  of  40  degrees  in  order  to  align  them  with 
the  strike  of  the  fault.   For  the  region  in  which  the 
tangential  components  had  not  been  rotated,  comparison  of  the 
Dixie  Valley  results  in  Table  8  showed  a  marked  agreement  in 
the  order  of  the  significance  of  the  results.   Therefore, 
these  interpretations  of  the  results  could  have  been  made 
from  a  single  analysis  of  the  rate  quantities  without  doing 
the  difference  calculations. 
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K-S   Statistic 
Component  N-S W-E 


Gravitational  Tide 
Normal  Tide 
Axial  Tide 
Tangent  Magnitude 
Total  Magnitude 

Gravitational  Tide  Rate 
Normal  Tide  Rate 
Axial  Tide  Rate 
Tangent  Magnitude  Rate 
Total  Magnitude  Rate 


.77 

.44 

.33 

.25 

.39 

.28 

.32 

.23 

.39 

.25 

42 

.36 

21 

.20 

25 

.23 

34 

.29 

53 

.42 

Results  of  Bias  Check  in  Tidal  Comparisons 
N-S  and  W-E  Extremes  of  the  San  Andreas  Area 
Using  the  K-S  Statistic 

Table  4 
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K-S  Statistic 
Component  N-S W-E 


Gravitational  Tide 
Normal  Tide 
Axial  Tide 
Tangent  Magnitude 
Total  Magnitude 


Gravitational  Tide  Rate  .46  .17 

Normal  Tide  Rate  .29  .17 

Axial  Tide  Rate  .26  .15 

Tangent  Magnitude  Rate  .30  .20 

Total  Magnitude  Rate  .61  .16 


.91 

.17 

.28 

.16 

.48 

.20 

.32 

.16 

.36 

.15 

Results  of  Bias  Check  in  Tidal  Comparisons 
N-S  and  W-E  Extremes  of  the  Dixie  Valley  Area 
Using  the  K-S  Statistic 

Table  5 
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Component  Statistic 

Gravitational  Tide  1.1411 

Normal  Tide  1.5038 

Axial  Tide  1.1776 

Tangent  Magnitude  1.3344 

Total  Magnitude  0.8658 

Gravitational  Tide  Rate  1.6515 

Normal  Tide  Rate  0.6126 

Axial  Tide  Rate  0.8469 

Tangent  Magnitude  Rate  0.9541 

Total  Magnitude  Rate  1.2273 

Gravitational  Tide  Difference  1.1724 

Normal  Tide  Difference  1.2752 

Axial  Tide  Difference  0.9117 

Tangent  Magnitude  Difference  1.2259 

Total  Magnitude  Difference  1.1985 


Results  of  Actual  Tidal  Comparisons  in  the 
San  Andreas  Area  Using  the  K-S  Statistic 

Table  6 
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Component 

Gravitational  Tide 
Normal  Tide 
Axial  Tide 
Tangent  Magnitude 
Total  Magnitude 


Statistic    Degrees  of  Freedom 


40.22 

19 

96.84 

18 

36.99 

17 

23.37 

18 

34.33 

18 

Gravitational  Tide  Rate 
Normal  Tide  Rate 
Axial  Tide  Rate 
Tangent  Magnitude  Rate 
Total  Magnitude  Rate 


114.21 
38.38 
23.48 
38.21 
81.64 


19 

18 
17 
18 
17 


Results  of  Actual  Tidal  Comparisons  in  the 
San  Andreas  Area  Using  the  Chi-Square  Statistic 

Table  7 
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Component  Statistic 

Gravitational  Tide  1.2218 

Normal  Tide  0.6278 

Axial  Tide  0.8792 

Tangent  Magnitude  0.7129 

Total  Magnitude  0.6139 

Gravitational  Tide  Rate  0.6535 

Normal  Tide  Rate  1.4119 

Axial  Tide  Rate  0.5901 

Tangent  Magnitude  Rate  0.4673 

Total  Magnitude  Rate  0.9228 

Gravitational  Tide  Difference  0.5743 

Normal  Tide  Difference  1.2654 

Axial  Tide  Difference  0.6040 

Tangent  Magnitude  Difference  0.5822 

Total  Magnitude  Difference  1.1605 


Results  of  Actual  Tidal  Comparisons  in  the 
Dixie  Valley  Area  Using  the  K-S  Statistic 

Table  8 
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Component 

Gravitational  Tide 
Normal  Tide 
Axial  Tide 
Tangent  Magnitude 
Total  Magnitude 


Statistic   Degrees  of  Freedom 


11.22 

13 

6.90 

14 

26.30 

14 

15.04 

15 

10.47 

11 

Gravitational  Tide  Rate 
Normal  Tide  Rate 
Axial  Tide  Rate 
Tangent  Magnitude  Rate 
Total  Magnitude  Rate 


7.34 

25.47 

5.99 

3.91 

13.72 


14 
14 
15 
12 
13 


Results  of  Actual  Tidal  Comparisons  in  the 
Dixie  Valley  Area  Using  the  Chi-Square  Statistic 

Table  9 
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VII.   CONCLUSIONS 

Results  of  the  analysis  indicate  that,  even  at  high  levels 
of  significance,  there  exist  components  of  the  Earth  tides 
and  their  rates -of -change  that  are  significantly  related  to 
the  occurrence  of  earthquakes  and  that  these  components  vary 
with  the  nature  and  orientation  of  the  fault  being  investi- 
gated. 

The  significance  of  the  normal  tidal  component  on  the 
strike-slip  San  Andreas  fault  gives  credence  to  the  theory 
that  the  occurrence  of  earthquakes  in  this  fault  zone  may 
result  from  the  lessening  of  the  frictional  force  acting 
between  the  two  vertical  faces  of  the  fault.   Since  the 
frictional  force  acting  between  two  surfaces  is  a  function  of 
the  normal  force  acting  on  those  surfaces,  a  lessening  of  the 
normal  force  could  decrease  the  frictional  force  enough  to 
allow  the  accumulated  stress  in  the  fault  to  be  released. 

The  relationship  between  the  rate-of -change  of  a  tidal 
force  and  the  difference  in  a  force  over  the  breadth  of  a 
fault  has  been  discussed  previously.   Although  no  difference 
quantities  were  considered  significant  at  a  0.05  level  of 
significance,  reference  to  Tables  2,  6  and  8  reveal  that  these 
quantities  did  differ  considerably  and  could  have  been  consid- 
ered significant  at  lower  levels  of  significance.   The  signif- 
icance of  the  gravitational  rate  in  the  region  of  the  San 
Andreas  fault  indicates  a  differential  in  the  vertical  forces 
acting  on  opposite  sides  of  the  fault.   This  significance  is 
difficult  to  interpret  in  view  of  the  strike-slip  nature  of 
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the  fault.   Because  of  the  proximity  of  this  fault  to  the 
Pacific  Ocean,  this  significance  may  indicate  the  presence  of 
ocean  tidal  loading  effects  on  the  occurrence  of  earthquakes 
in  this  region.   The  significance  of  the  normal  rate  in  the 
region  of  the  Dixie  Valley  fault  indicates  a  differential  in 
the  east-west  force  acting  on  opposite  sides  of  this  fault. 
This  interpretation  is  consistent  with  the  nature  of  the 
Dixie  Valley  fault  in  that  the  downward  slippage  of  the  center 
of  the  fault  zone  occurs  with  the  separation  of  the  sides 
of  the  fault. 
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APPENDIX  A 


A.  OBSERVATIONS  ON  ASSOCIATED  LEGENDRE  POLYNOMIALS. 

The  generating  function  for  the  associated  Legendre 

th 
polynomial  of  the  m   order  is 

Other  useful  relationships  are 


X 


(o 


P^M-O-*)*^  P.W  (*> 


R+.W--  TTrfcz+OMPiM-l-Pt-X*)]  O) 


Letting  m=l  in  (1)  and  rearranging  terms, 

3  1    w 

(.-»**  *yl)~*  =  (i-^'^v'"-^  ^ 

From    (2) , 

(,-**)- '/l  P.'Cy)  r    i  o> 

0-*^~*  P*'M  =  3x 

« 

From    (3) , 

P.  00*  I  <y 

P,  C*)  =  x 

P,Wri(3x'-0 

P4  C*)  =  •g'C35x''-3oxZ-e3) 


32 


B. 


EARTH  TIDE  CALCULATIONS 


R  =  Earth  radius  vector  to  epicenter 

a  =  Scalar  magnitude  of  R 

R-   =  Distance  vector  from  epicenter  to  B- 

r-   =  Scalar  magnitude  of  R- 

R-   =  Distance  vector  from  origin  to  B- 

r-   =  Scalar  magnitude  of  R- 

B-  =  Extra-terrestrial  body,  i-1  (moon),  2  (Sun) 

F-   =  Force  vector  between  the  origin  and  B- 

F-   =  Force  vector  between  the  origin  and  B- 
—10  &        l 

tf'i  =  Angle  between  R  and  R- 
m-  =  mass  of  B- 


From  Newton's  laws  of  gravitation, 


Fie- 


F 


CO 


K 


The  total  tidal  force  vector,  T,  is  the  sum  over  all  bodies 


T  -  .1  Ti 

where  T^  =  F^0  -  F_^e  =  k*m^ 


*io   he 


3  3 

r .     r  - 
io    le 


Since 


k-  r*u 


V-c« 


With  l^e  =  ho  +<*.  -Za«ir;0  c©s  (<**.),   the  expressi 
be  written  as 


on 


V-L, 


Ke 


to 


can 
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h« 


Ke, 


^2(^o)C°S(k:) 


q 


J3 


(?) 


Letting   x^cos  (c*;. )  ,    and  yi  =   -§-   >    a   comparison   of    (4)    with 


(8)    yields, 
v3 


10 

2    -—         »•'      ,-.  / 


00 


?-'_i R(c-5<0) 


oo 


c(coS(oO 


Substituting    (9)    into    (7)    yields, 

~h7~ 


Ti 


.  „s~  /•-5_A5"'__i__R(c0s(0<;)) 


Letting  t   be  the  component  of  T^  in  the  direction  of  N,  a 
unit  vector  in  the  plane  normal  to  R,  and  letting  (3.  be  the 
angle  between  N  and  T • ,  results  in  the  general  expression, 
tn=N-T- ,  and  after  rearranging  terms, 


k*  i^;. 


?■'  j 


Expanding  terms  in  (10)  and  substituting  values  from  (5) , 


K  •  v^i 


Z~   COS 


io 


0»;) 


fef 


d 


Noting 'that    cos(o<;)  = 


R-R 


10 


xi 


a*  r •  r • 

10  10 


if   N    is    taken    to   be   the   north(z)    direction,    (11)    gives    t, 
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w 


ith  cos(ft)= 


z- 


10 


,  and  if  N  is  taken  to  be  the  east(y) 


direction,  (11)  gives  t   with  cos(/^l)  = 


y- 


y 


10 


If  the 


co-ordinate  system  is  oriented  so  that  R  is  in  the  x-direc 

R-T. 
tion,  then  tv=t  = 


— 1 


CO 


+  1    fe) 


~  co3Co<0 


Expanding  a  few  terms  of  (12)  ,  collecting  coefficients  of 


a  U 


where  j  =0 , 1 ,  and  substituting  values  from  (6), 
3c©560-(  +f-ryr-)[S"cos  (o(i)-3cos(o< 


with  cos(p<i)  = 


co  X,  u 


;>H 


10 


For  tangential  components  in  other  than 


the  y  and  z  directions ,  N  can  be  any  unit  vector  in  the 
plane  normal  to  the  vector  R. 
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